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ABSTRACT Male circumcision reduces female-to-male HIV transmission. Hypothesized mechanisms for this protective effect 
include decreased HIV target cell recruitment and activation due to changes in the penis microbiome. We compared the coronal 
sulcus microbiota of men from a group of uncircumcised controls (n = 77) and from a circumcised intervention group (n = 79) 
at enrollment and year 1 follow-up in a randomized circumcision trial in Rakai, Uganda. We characterized microbiota usingl6S 
rRNA gene-based quantitative PCR (qPCR) and pyrosequencing, log response ratio (LRR), Bayesian classification, nonmetric 
multidimensional scaling (nMDS), and permutational multivariate analysis of variance (PerMANOVA). At baseline, men in 
both study arms had comparable coronal sulcus microbiota; however, by year 1 , circumcision decreased the total bacterial load 
and reduced microbiota biodiversity. Specifically, the prevalence and absolute abundance of 12 anaerobic bacterial taxa de- 
creased significantly in the circumcised men. While aerobic bacterial taxa also increased postcircumcision, these gains were mi- 
nor. The reduction in anaerobes may partly account for the effects of circumcision on reduced HIV acquisition. 

IMPORTANCE The bacterial changes identified in this study may play an important role in the HIV risk reduction conferred by 
male circumcision. Decreasing the load of specific anaerobes could reduce HIV target cell recruitment to the foreskin. Under- 
standing the mechanisms that underlie the benefits of male circumcision could help to identify new intervention strategies for 
decreasing HIV transmission, applicable to populations with high HIV prevalence where male circumcision is culturally less ac- 
ceptable. 
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Male circumcision (MC) reduces the risk of HIV acquisition in 
men by 50 to 60% (1-3) and decreases the incidence and 
prevalence of herpes simplex virus 2 (HSV-2) (4) and human pap- 
illomavirus (HPV) (4, 5). The impact of MC on classical bacterial 
sexually transmitted infections (STIs), such as Chlamydia tracho- 
matis, Neisseria gonorrhoeae, Treponema pallidum, and Trichomo- 
nas vaginalis infection, is more equivocal (4, 6-8). Women with 
circumcised male partners are at lower risk for STIs ranging from 
HPV to Trichomonas vaginalis infection (6, 9). This suggests that 
MC reduces the risk of viral STIs in men and of STI transmission 
to their female partners (10). 

MC is hypothesized to reduce HIV risk in men by changing the 
penile anatomy and by altering the genital microbiology (11). 
With respect to the anatomic changes, MC removes the prepuce, 
which decreases the number of available HIV target cells on the 
penis (11, 12). It remains unclear whether decreases in viral STIs 
post-MC contribute to HIV risk reduction. HSV-2 infection in- 



creases the risk of HIV in observational studies (13, 14), but trials 
aimed at controlling viral and classical bacterial STIs have largely 
failed to reduce HIV transmission (15, 16). Removal of the prepu- 
tial tissue also eliminates the moist subpreputial environment, 
which can modify the genital bacterial communities (i.e., the mi- 
crobiota) and may have a broad impact on the genital microbiol- 
ogy (17). 

Recently, genital epithelial inflammation associated with bac- 
terial antigens has emerged as a possible factor in increasing sus- 
ceptibility of genital HIV target cells (18-23). These findings sug- 
gest that specific groups of genital bacteria, including those not 
associated with classical STIs, could elicit local immune responses 
that promote epithelial inflammation and recruitment of HIV tar- 
get cells. Thus, changes in the genital bacterial microbiota could be 
linked to HIV acquisition. 

Previously, we reported the impact of MC on the coronal sul- 
cus microbiota composition in 12 men (17). However, this study 
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TABLE 1 Demographic characteristics, sexual behaviors, and symptoms of sexually transmitted infections for the control and intervention arms at 
enrollment 

No. (%) in group 

Characteristic Control (n = 77) Intervention (n = 79) 



Age (yr) 



15-19 


1 fl 3) 


1 (1 3) 


20-24 


13 (16.9) 


15 (19.0) 


25-29 


21 (27.3) 


25 (31.6) 


30^49 


42 (54.5) 


38 (48.1) 


A'larital status 






Currently married, monogamous 


70 (91.0) 


71 (89.9) 


Currently married, polygamous 


7 (9.0) 


8 (10.1) 


No. of sexual partners in past yr 






1 


45 (58.4) 


45 (57.0) 


2 


23 (29.9) 


26 (32.9) 


>3 


9 fll 7) 


8 (10.1) 


Nonmarital sexual relationships 






No 


66 (85.7) 


67 (84.8) 


Yes 


11 (14.3) 


12 (15.2) 


Condom use in past yr 








46 (59.7) 


54 (68.3) 


Inconsistent use 


30 (39.0) 


24 (30.4) 


Consistent use 


1(1.3) 


1(1.3) 


Syphilis infection 






No 


72 (93.5) 


73 (92.4) 


Yes 


4 (5.2) 


5 (6.3) 


Not tested 


1(1.3) 


1(1.3) 


HSV-2 infection 






No 


36 (46.8) 


37 (46.8) 


Yes 


32 (41.5) 


32 (40.5) 


Indeterminate 


9(11.5) 


9(11.4) 


Not tested 


1(1.3) 


1(1.3) 


Self-reported symptoms of sexually transmitted infection 






Genital ulcer disease 


0 (0.0) 


0 (0.0) 


Urethral discharge 


0 (0.0) 


1(1.3) 


Dysuria 


0 (0.0) 


0 (0.0) 



lacked uncircumcised controls. In the current study, we assessed 
the effect of MC on the genital microbiota using absolute abun- 
dance. In addition, we applied novel analyses to assess the micro- 
biota changes attributable to MC. We hypothesized that MC 
would significantly decrease coronal sulcus bacterial abundance 
and modify the microbiota in participants randomly assigned to 
receive MC but not in those who remained uncircumcised. Here, 
we report a study of penile coronal sulcus microbiota in 77 control 
and 79 intervention-arm participants from the Rakai MC ran- 
domized controlled trial in Uganda. 

RESULTS 

Study participant profile at enrollment. At enrollment, men 
from the control and intervention arms had similar sociodemo- 
graphic characteristics, sexual practices, sexually transmitted in- 
fections, and symptoms (Table 1). 

Coronal sulcus bacteria in the uncircumcised penis at enroll- 
ment, (i) Prevalence. At enrollment, the prevalences of coronal 
sulcus bacterial were comparable between the two study arms (Ta- 
ble 2). Some of the most common coronal sulcus bacteria seen at 
enrollment included those from the Prevotellaceae, Veillonellaceae, 
Clostridiales family XI, Actinomycetaceae, Coriobacteriaceae, and 
Porphyromonadaceae. Two groups of bacteria from the order Clos- 
tridiales were highly prevalent but could not be assigned with suf- 
ficient confidence to known lower taxa and are referred to as un- 



classified Clostridiales family XI and unclassified Clostridiales 
(Table 2). 

(ii) Relative abundance. Most coronal sulcus bacteria were 
observed in relatively low abundances (Table 2). Prevotella spp. 
were the most dominant, followed by unclassified members of the 
Clostridiales and Cory neb acterium spp. Six others — Peptoniphilus 
spp., Anaerococcus spp., Finegoldia spp., Murdochiella spp., Por- 
phyromonas spp., and Lactobacillus spp. — were found at relative 
abundances of approximately 5%. The remaining coronal sulcus 
bacteria were detected at lower than 1% (Table 1). 

Male circumcision reduces coronal sulcus bacterial load. At 
enrollment, similar mean bacterial loads were seen in the two 
study groups based on measurements of the bacterial 16S rRNA 
gene, with an average of 1.4 X 10 5 copies (standard deviation [SD] 
= 3.1 X 10 5 ) in the control arm and 2.0 X 10 5 copies (SD = 4.8 X 
10 5 ) in the intervention arm. At year 1, the total bacterial load 
decreased significantly in both arms. In the uncircumcised men, 
the average bacterial load decreased to 5.7 X 10 4 copies (SD = 1.19 
X 10 5 ), but the circumcised men had an average of 3.8 X 10 4 
copies (SD = 1.80 X 10 5 ) (log response ratio P = 0.048) (Fig. 1). 
Thus, MC significantly decreased the coronal sulcus bacterial load 
relative to changes in uncircumcised men. 

Male circumcision significantly altered prevalences of coro- 
nal sulcus bacteria. Fifteen coronal sulcus bacteria significantly 
decreased in prevalence post-MC (P < 0.05), among which 12 
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TABLE 2 Prevalences and proportional abundances of the 40 most common coronal sulcus bacteria in the control and intervention arms at 
enrollment" 



Family 


Genus 


Prevalence (%) in 


group 


Avg proportional abundance [% (SD)] in group 


Control (« = 77) 


Intervention (n = 79) 


Control 


Intervention 


Clostridiales family XIa 


Peptoniphilus spp. 


74 (96.1) 


72 (91.1) 


5.4 (5.9) 


5.1 (5.8) 


Clostridiales family XIa 


Anaerococcus spp. 


71 (92.2) 


68 (86.1) 


5.1 (8.1) 


4.3 (6.2) 


NA 


Unclassified Clostridiales 


69 (89.6) 


70 (88.6) 


15.9 (16.0) 


14.3 (14.5) 


Prevotellaceae 


Prevotella spp. 


69 (89.6) 


67 (84.8) 


21.4 (17.0) 


23.1 (20.5) 


Clostridiales family XIa 


Finegoldia spp. 


63 (81.8) 


64(81.0) 


6.5 (8.3) 


7.1 (10.6) 


Clostridiales family XIa 


Murdochiella spp. 


62 (80.5) 


58 (73.4) 


3.2 (4.7) 


5.2 (8.6) 


Porphyromonadaceae 


Porphyromonas spp. 


61 (79.2) 


55 (69.6) 


5.4 (6.1) 


4.8 (8.4) 


Corynebacteriaceae 


Corynebacterium spp. 


57 (74.0) 


52 (65.8) 


12.2 (21.2) 


8.5 (18.8) 


Clostridiales family XI 


Unclassified Clostridiales Family XI 


54 (70.1) 


53 (67.1) 


0.7 (0.8) 


0.8 (1.5) 


Veillonellaceae 


Dialister spp. 


53 (68.8) 


43 (54.4) 


1.6 (1.9) 


1.0 (1.6) 


Veillonellaceae 


Negativicoccus spp. 


40 (51.9) 


36 (45.6) 


1.0 (2.4) 


0.8 (1.8) 


Peptostreptococcaceae 


Peptostreptococcus spp. 


31 (40.3) 


36 (45.6) 


0.9 (2.1) 


1.3 (2.7) 


Actinomycetaceae 


Mobiluncus spp. 


38 (49.4) 


26 (32.9) 


1.5 (5.2) 


1.1 (3.8) 


Bifidobacteriaceae 


Gardnerella spp. 


33 (42.9) 


25 (31.6) 


1.8 (5.9) 


0.9 (3.2) 


Lactobacillaceae 


Lactobacillus spp. 


26 (33.8) 


28 (35.4) 


2.4 (10.7) 


8 (21.0) 


Staphylococcaceae 


Staphylococcus spp. 


29 (37.7) 


21 (26.6) 


1.4 (3.9) 


1.1 (6.0) 


Ruminococcaceae 


Saccharofermentans spp. 


28 (36.4) 


21 (26.6) 


0.3 (0.6) 


0.5 (1.9) 


Streptoeoccaceae 


Streptococcus spp. 


26 (33.8) 


19 (24.1) 


1.3 (6.1) 


0.5 (2.8) 


Actinomycetaceae 


Actinomyces spp. 


26 (33.8) 


17 (21.5) 


0.1 (0.5) 


0.1 (0.2) 


Veillonellaceae 


Veillonella spp. 


26 (33.8) 


14 (17.7) 


0.5 (1.5) 


0.6 (3.1) 


Peptococcaceae 1 


Peptococcus spp. 


25 (32.5) 


14 (17.7) 


0.1 (0.1) 


0.04 (0.1) 


Coriobacteriaceae 


Olsenella spp. 


20 (26.0) 


19 (24.1) 


0.1 (0.2) 


0.1 (0.2) 


Actinomycetaceae 


Arcanobacterium spp. 


23 (29.9) 


14 (17.7) 


0.1 (0.2) 


0.1 (0.2) 


Lachnospiraceae 


Howardella spp. 


19 (24.7) 


9 (11.4) 


0.1 (0.1) 


0.03 (0.1) 


Clostridiales family XIa 


Parvimonas spp. 


17 (22.1) 


10 (12.7) 


0.2 (0.5) 


0.3 (1.2) 


Coriobacteriaceae 


Atopobium spp. 


14 (18.2) 


12 (15.2) 


0.1 (0.2) 


0.2 (1.1) 


Leptotrichiaceae 


Sneathia spp. 


13 (16.9) 


13 (16.5) 


0.2 (0.8) 


0.3 (1.4) 


Sutterellaceae 


Sutterella spp. 


13 (16.9) 


12 (15.2) 


0.1 (0.2) 


0.03 (0.1) 


Lachnospiraceae 


Moryella spp. 


14 (18.2) 


7 (8.9) 


0.1 (0.2) 


0.03 (0.1) 


Peptostreptococcaceae 


Peptostreptococcaceae family 


12 (15.6) 


9 (11.4) 


0.1 (0.5) 


0 03 (0 11 


Spirochaetaceae 


Treponema spp. 


10 (13.0) 


11 (13.9) 


0.2 (0.6) 


0.2 (0.5) 


Fusobacteriaceae 


Fusobacterium spp. 


8 (10.4) 


13 (16.5) 


0.2 (0.8) 


0.8 (3.8) 


Synergistaceae 


Pyramidobacter spp. 


13 (16.9) 


7 (8.9) 


0.2 (0.8) 


0.2 (0.8) 


Aerococcaceae 


Facklamia spp. 


12 (15.6) 


8 (10.1) 


0.1 (0.3) 


0.1 (0.5) 


Clostridiales family XIa 


Anaerosphaera spp. 


9 (11.7) 


11 (13.9) 


0.02 (0.1) 


0.1 (0.2) 


Micrococcaceae 


Kocuria spp. 


10 (13.0) 


8 (10.1) 


0.05 (0.2) 


0.1 (0.2) 


Veillonellaceae 


Megasphaera spp. 


10 (13.0) 


8 (10.1) 


0.2 (0.6) 


0.1 (0.4) 


Micrococcaceae 


Micrococcus spp. 


8 (10.4) 


10 (12.7) 


0.04 (0.2) 


0.04 (0.1) 


Bacillales family XI 


Gemella spp. 


10 (13.0) 


6 (7.6) 


0.04 (0.1) 


0.03 (0.1) 


Burkholderiaceae 


Ralstonia spp. 


13 (16.9) 


2 (2.5) 


0.1 (0.3) 


0.01 (0.1) 



a *, false discovery rate (FDR)-adjusted P value < 0.05. 



were strict anaerobes, including Porphyromonas spp. 
(AAPrevalence = —43.10%), Prevotella spp. (AAPrevalence = 
-34.21%), Negativicoccus spp. (AAPrevalence = -28.95%), Di- 
alister spp. (AAPrevalence = —30.18%), Mobiluncus spp. 
(AAPrevalence = —13.69%), and six genera from Clostridiales 
family XI, among others (Table 3). The reductions in anaerobe 
prevalence due to MC were often substantial, but MC did not 
significantly reduce all anaerobes; notably, Atopobium spp., 
Sneathia spp., and Megasphaera spp. showed no statistically sig- 
nificant decrease post-MC. 

Seven coronal sulcus bacteria became more prevalent post- 
MC. Among these, five also increased in prevalence in the uncir- 
cumcised men over time, suggesting either an effect of time or 
changes in behavior with participation in the trial. Nevertheless, a 
greater number of the circumcised than of the uncircumcised men 
acquired these specific bacteria, as shown by the positive 
AAPrevalence values (Table 3). The aerobic Kocuria spp. and the 



facultative anaerobic Facklamia spp. were the two types of bacteria 
that became more prevalent exclusively in the circumcised men. 
Other bacteria were uncommon in the uncircumcised penis but 
increased in prevalence post-MC (see Table SI in the supplemen- 
tal material). 

Male circumcision modified coronal sulcus microbiota bio- 
diversity and composition, (i) Microbiota biodiversity. MC sig- 
nificantly reduced the evenness of the microbiota, indicating a 
general decrease in the number of dominant coronal sulcus bac- 
teria post-MC (E treatment effect = -0.053; 95% CI = -0.101 to 
—0.005). In addition, MC also significantly decreased the struc- 
tural diversity of the microbiota (D treatment effect = —1.26; 
95% CI = -2.04 to -0.52) 

(ii) Microbiota composition. MC reshaped the composition 
of the coronal sulcus microbiota, producing a more homogeneous 
post-MC profile (Fig. 2A and B). While both study groups show 
significant temporal changes in microbiota composition, the 
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FIG 1 Changes in the coronal sulcus bacterial load as measured by the log 
response ratio for the uncircumcised (Contr.; in red) versus the circumcised 
(Interv.; in orange) men, shown by group (top panel) and by individual (bot- 

(Continued) 



change was more marked in the circumcised men (PerMANOVA 
F statistic = 13.1; P < 0.001) (Fig. 2A; see also Fig. SI in the 
supplemental material) than in the uncircumcised men (Per- 
MANOVA F statistic = 3.11 P = 0.02) (Fig. 2B; see also Fig. SI). 

Circumcision significantly reduced previously abundant 
coronal sulcus bacteria. To quantify the impact of MC on coronal 
sulcus bacteria, we determined the MC effect size. This was per- 
formed for genera that either significantly decreased (i.e., "nega- 
tive responders") or increased (i.e., "positive responders") after 
MC (Table 4). Among the negative responders, Prevotella spp., 
Porphyromonas spp., Finegoldia spp., and Peptostreptococcus spp. 
decreased in both prevalence and absolute abundance, with effec- 
tive load reductions ranging from — 1,157 to —25,327 16S rRNA 
gene copies (Table 4). Other negative responders decreased signif- 
icantly in either prevalence [n = 8) or absolute abundance (n = 2). 
Several negative responders had substantial effective load reduc- 
tions that were also highly variable, such as unclassified Clostridi- 
als, Peptoniphilus spp., and Murdochiella spp. As a result, they 
have large but non-statistically significant effect sizes (Table 4). 

In contrast, many positive responders had smaller MC effect 
sizes that were statistically significant, a finding that indicated a 
more uniform bacterial gain among circumcised men. On aver- 
age, Corynebacterium spp. increased by 2,860 and Staphylococcus 
spp. by 249 16S rRNA gene copies per individual (Table 4). The 
third-highest mean MC effect size was seen in Helcococcus spp., 
which belong to Clostridiales family XI, and that response con- 
trasts with the broadly negative impact of MC on other Clostridi- 
ales family XI members. Overall, the relatively larger MC effect 
sizes in negative responders indicate that MC primarily reduced 
previously abundant coronal sulcus bacteria, accompanied by 
other minor abundance gains. 

DISCUSSION 

In a randomized trial of MC, we showed that MC significantly 
reduced the bacterial load by reducing both the prevalence and the 
absolute abundance of many coronal sulcus bacteria. The two 
study groups had comparable coronal sulcus microbiota at enroll- 
ment that consisted of multiple microbiota types, but MC pro- 
foundly altered the composition of the microbiota and reduced its 
biodiversity. Over time, changes in the coronal sulcus microbiota 
were observed in the uncircumcised men. However, after adjust- 
ing for these temporal changes, we found that there were signifi- 
cantly greater decreases in the total bacterial load, microbiota bio- 
diversity, and microbiota composition in the circumcised men 
that were attributable to MC. 

The role of coronal sulcus bacteria in heterosexual HIV acqui- 
sition remains unknown. Recent studies suggest that the non-STI 
genital bacteria may affect the susceptibility of foreskin HIV target 
cells (22, 24). Of the HIV target cell types found in the foreskin, 



Figure Legend Continued 

torn panel). In the group comparison, the box of each box plot denotes the 
interquartile range (IQR) (quartile 1 [Ql] to Q3) and the corresponding me- 
dian, whereas the whiskers signify the upper and lower 1 .5 X IQR. Outliers are 
shown as open symbols in each box plot. There was a statistically significant 
reduction in bacterial load for the circumcised men compared to that for the 
uncircumcised men (P = 0.048). As shown in the scatter plot in the bottom 
panel, although a decrease was observed for many individuals from both 
groups, more circumcised men showed decreases (i.e., negative log response 
ratios) (62/79, 78.5%) than did those that remained uncircumcised (51/77, 
66.2%). 
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TABLE 3 Prevalences and changes in 


prevalence of the 40 most common 


coronal sulcus bacteria for 


uncircumcised and circumcised men at year 1 






Prevalence (yo) 6 






Bacterial group 


Oxygen tolerance* 1 


TT -j 

Uncircumcised 


Circumcised 


AAPrevalence (%)' 


Peptoniphilus spp. 


AN 


71 (92.2) 


37 (46.8)** 


-40.41 


Anaerococcus spp. 


AN 


71 (92.2) 


59 (74.7) 


-11.39 


Unclassified Clostridiales spp. 


NA 


65 (84.4) 


38 (48.1)** 


-35.31 


Prevotella spp. 


AN 


70 (90.9) 


41 (51.9)** 


-34.21 


Finegoldia spp. 


AN 


61 (79.2) 


38 (48.1)** 


-30.31 


Murdochiella spp. 


AN 


51 (66.2) 


12 (15.2)** 


-43.94 


Porphyromonas spp. 


AN 


62 (80.5) 


23 (29.1)** 


-43.10 


Corynebacterium spp. 


FAN 


71 (92.2) 


77 (97.5)** 


13.46 


Unclassified Clostridiales family XI 


NA 


55 (71.4) 


20 (25.3)** 


-43.07 


Dialister spp. 


AN 


47 (61.0) 


13 (16.5)** 


-30.18 


Negativicoccus spp. 


AN 


35 (45.5) 


8 (10.1)** 


-28.95 


Peptostreptococcus spp. 


AN 


32 (41.6) 


13 (16.5)* 


-30.41 


Mobiluncus spp. 


AN 


31 (40.3) 


8 (10.1)* 


-13.69 


Gardnerella spp. 


FAN 


34 (44.2) 


32 (40.5) 


7.56 


Lactobacillus spp. 


FAN/AN/MAE 


32 (41.6) 


35 (44.3) 


1.07 


Staphylococcus spp. 


FAN 


48 (62.3)* 


69 (87.3)** 


36.08 


Saccharofermentans spp. 


AN 


21 (27.3) 


10 (12.7) 


-4.83 


Streptococcus spp. 


FAN 


34 (44.2) 


28 (35.4) 


1.00 


Actinomyces spp. 


FAN 


27 (35.1) 


8 (10.1) 


-12.69 


Veillonella spp. 


AN 


26 (33.8) 


16 (20.3) 


2.53 


Peptococcus spp. 


AN 


14 (18.2) 


7 (8.9) 


5.42 


Olsenella spp. 


AN 


19 (24.7) 


5 (6.3) 


-16.42 


Arcanobacterium spp. 


FAN 


15 (19.5) 


4(5.1) 


-2.27 


Howardella spp. 


AN 


18 (23.4) 


5 (6.3) 


-3.76 


Parvimonas spp. 


AN 


13 (16.9) 


12 (15.2) 


7.73 


Atopobium spp. 


AN 


18 (23.4) 


11 (13.9) 


-6.46 


Sneathia spp. 


AN 


15 (19.5) 


9 (11.4) 


-7.66 


Sutterella spp. 


AN 


15 (19.5) 


2 (2.5)* 


15.26 


Moryella spp. 


AN 


11 (14.3) 


4 (5.1) 


0.10 


Peptostreptococcaceae family 


NA 


4 (5.2) 


1 (1.3) 


0.26 


Treponema spp. 


AN 


8 (10.4) 


4(5.1) 


-6.26 


Fusobacterium spp. 


AN 


17 (22.1) 


11 (13.9) 


-14.22 


Pyramidobacter spp. 


AN 


7 (9.1) 


3 (3.8) 


2.73 


Facklamia spp. 


FAN 


17 (22.1) 


21 (26.6) 


9.96 


Anaerosphaera spp. 


AN 


9 (11.7) 


7 (8.9) 


-5.06 


Kocuria spp. 


AE 


20 (26.0) 


36 (45.6)** 


22.46 


Megasphaera spp. 


AN 


11 (14.3) 


10 (12.7) 


1.23 


Micrococcus spp. 


AE 


22 (28.6) 


39.24 


8.40 


Gemella spp. 


FAN 


15 (19.5) 


13 (16.5) 


2.37 


Ralstonia spp. 


AE 


16 (20.8) 


2 (2.5) 


-3.90 



a AN, strictly anaerobic; AE strictly aerobic; FAN, facultative anaerobic; MAE, microaerophilic; NA, no data. 

b **, FDR-adjusted P value < 0.0001 for change in prevalence over time (i.e., APrevalence); *, FDR-adjusted P value < 0.05 for APrevalence. 

c AAPrevalence, shown as a percentile, is the change in prevalence seen for the circumcised men over time compared to that for men that remain uncircumcised. 



Langerhans cells (LCs) have been hypothesized to play a key role 
in mediating HIV infection (24). Located proximal to the epithe- 
lial surface, naive LCs bind, internalize, and degrade HIV parti- 
cles; however, when activated by a high HIV load, active STIs, or 
bacterium-associated inflammatory mediators, such as lipopoly- 
saccharide (LPS) and tumor necrosis factor alpha (TNF-a), LCs 
bind and present HIV particles to CD4+ T cells (24, 25). 

As already mentioned, the changes in penile microbiota and 
STI incidence after MC may be attributable not only to the ana- 
tomic alteration itself but also to behavioral changes in circum- 
cised men or men enrolled in a clinical trial in general. However, 
analysis of the Rakai data showed that MC did not significantly 
alter behavior during the trial ( 1 ) . Likewise, the natural dynamics 
of the circumcised coronal sulcus microbiota are unknown, but 
sampling performed 1 year postcircumcision likely represents a 
persistent change. Our use of novel analysis metrics, such as the 
log response ratio and MC effect size, permitted adjustment for 



the impact of time and trial participation, which allowed us to 
quantify microbiota changes attributable solely to MC. 

MC has been associated with reduction of bacterial vaginosis 
(BV) in female sexual partners, but the sharing of genital micro- 
biota between sexual partners is not well understood (6, 26). We 
found that a subset of bacteria associated with BV decreased after 
MC, including Prevotella spp., Fusobacterium spp., and Mobilun- 
cus spp., while others, such as Gardnerella spp., Sneathia spp., 
Actinomyces spp., Atopobium spp., Megasphaera spp., and Veillo- 
nella spp., were not significantly altered. 

MC selected for bacteria capable of surviving in the aerated 
circumcised microenvironment. At enrollment, the microbiota 
types were comparable in the two study groups. However, nearly 
all of the bacteria that decreased after MC were strict anaerobes, 
except for Actinomyces spp. and Arcanobacterium spp., which are 
facultative anaerobes. We also show that as a facultative anaerobe 
(27), Helcococcus spp. constituted the single positive responder to 
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FIG 2 The nonmetric multidimensional scaling (nMDS) ordination plots 
enable the visualization of individuals' microbiota over time. In nMDS plots, 

(Continued) 



MC in Clostridiales family XI. The large competitive advantage of 
a single genus from a previously diverse and abundant bacterial 
family illustrates the strong and functionally cohesive selective 
pressure exerted by MC through changes to the coronal sulcus 
microenvironment. 

One of the largest positive responders to MC was the Staphy- 
lococcus species group. Although we did not attempt to perform 
species-level analysis, the most abundant Staphylococcus bacteria 
on the post-MC coronal sulcus included S. haemolyticus, S. homi- 
nis, S. epidermidis, S. xylosus, and their genetic near neighbors. It is 
important to note that S. aureus and S. epidermidis are common 
commensals on exposed human epithelial and mucosal surfaces. 
Thus, their increase post-MC is unlikely to affect the pathogenic 
potential. 

We integrated culture-independent bacterial identification, an 
ecological analytical framework, and a randomized study design 
to reveal the impact of MC on the penis microbiome. Combining 
bacterial quantification with parallel sequencing showed that cir- 
cumcision resulted in significant decreases in the absolute abun- 
dances of several anaerobic bacterial taxa that defined the uncir- 
cumcised penis microbiome. Currently, we know little about the 
role of these fastidious anaerobes in the male urogenital tract or 
the broader context of human health. Future studies are required 
to determine if a decreased anaerobic bacterial load modifies fore- 
skin inflammation and HIV target cell recruitment/susceptibility, 
which may play a role in HIV risk reduction conferred by MC. 

MATERIALS AND METHODS 

Study design and subjects. We conducted a randomized trial of MC for 
HIV prevention in 2004 to 2006 (1). In this study, HIV-negative, uncir- 
cumcised men of ages 15 to 49 were randomized to either immediate 
circumcision (intervention group) or circumcision delayed for 24 months 
(control group), as described previously (1, 4). All circumcision proce- 
dures were performed in one surgical facility by the same team of urolo- 
gists and trained medical officers, using a single surgical procedure, the 
"sleeve method" of circumcision. Prophylactic antibiotics were not used 
for these procedures, and antibiotic use in this geographic region was 
minimal. Study participants were provided access to regular reproductive 
health services and followed at 6, 12, and 24 months to assess HIV and 
sexually transmitted infection (STI) acquisition, as described in detail 
elsewhere ( 1 ) . Men with a diagnosis of syphilis or symptoms suggestive of 
an STI were treated, but this was uncommon. Specifically, there were three 
new cases of syphilis in the control group and one new case in the inter- 
vention group who were treated with intramuscular benzathine penicillin. 
In addition, five control group men but none from the intervention group 
were treated for STI symptoms at month 6, and three control group men 
and three intervention group men were treated at year 1 . Symptomatic 
patients were given azithromycin, ciprofloxacin, or metronidazole. 
HSV-2 and HPV data were not available at the time of the trial, and thus 
no treatment was given. 

At each visit, clinicians collected penile swabs from the coronal sulcus 
as follows. Sterile cotton-tipped applicators (Thermo, Fisher Scientific, 
Waltham, MA) were premoistened with sterile saline and rolled over the 
coronal sulcus twice in a nontraumatic fashion. The swabs were immedi- 
ately placed in 1 ml of Amplicor specimen transport medium (Roche 



Figure Legend Continued 

each data point represents an individual's microbiota at one time point. The 
centroids and 95% confidence ellipses for each study group are as shown. Here, 
the coronal sulcus microbiota in men that remained uncircumcised showed 
minor variations from enrollment (in blue) to year-1 (in orange) (Fig. 2A). In 
contrast, significant shifts were seen in the circumcised men (Fig. 2B). 
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TABLE 4 Effect size of MC, measured as the change in absolute abundances of coronal sulcus bacteria that significantly decreased (i.e., "negative 
responders") or increased (i.e., "positive responders") post-MC, adjusted by changes in abundance among uncircumcised men over time" 

Category and bacterial group 6 Indicator value (FDR-adjusted P value) MC effect size [mean (90% CI)] 

Negative responders to MC 



Prevotella spp.*** 


0.18 (0.06) 


-25,327 (-48,812 to -3,988) 


Porphyromonas spp.*** 


0.27 (<0.01) 


-14,232 (-28,698 to -2,358) 


Unclassified Clostridials spp.* 


0.24 (0.01) 


-10,087 (-32,278 to 12,802) 


Unclassified Clostridiales family XI* 


0.22 (0.03) 


-3,299 (-8,715 to 647) 


Murdochiella spp.* 


0.30 (<0.01) 


-3,207 (-10,314 to 5,319) 


Peptoniphilus spp.* 


0.29 (<0.01) 


-1,349 (-8,121 to 5,388) 


Finegoldia spp.*** 


0.24 (0.07) 


-1,343 (-2,385 to -438) 


Anaerococcus spp.** 


0.27 (0.03) 


-1,284 (-2,292 to -337) 


Peptostreptococcus spp.*** 


0.21 (<0.01) 


-1,157 (-2,415 to -188) 


Mobiluncus spp.* 


0.16 (0.02) 


-568 (-1,480 to 290) 


Actinomyces spp. 


0.20 (0.01) 


-286 (-780 to 13) 


Saccharofermentans spp. 


0.16 (0.03) 


-281 (-1,545 to 1,032) 


Negativicoccus spp.* 


0.28 (<0.01) 


-244 (-785 to 321) 


Sutterella spp. * 


0.13 (0.01) 


-110 (-258 to 12) 


Howardella spp. 


0.14 (0.01) 


-44 (-175 to 75) 


Olsenella spp.* 


0.17 (0.01) 


-37 (-270 to 171) 


Ralstonia spp. 


0.14 (0.01) 


0 (0 to 1) 


Dialister spp.* 


0.24 (<0.01) 


1,434 (-1,259 to 4,618) 


)sitive responders to MC 






Rothia spp. 


0.12 (<0.01) 


1 (0 to 2) 


Roseomonas spp. 


0.11 (0.04) 


1 (0 to 2) 


Kocuria spp.*** 


0.32 (<0.01) 


8 (5 to 11) 


Brevibacterium*** 


0.30 (<0.01) 


13 (5 to 22) 


Eremococcus** 


0.26 (0.06) 


45 (8 to 94) 


Bifidobacterium 


0.16 (<0.01) 


101 (-7 to 243) 


Helcococcus*** 


0.21 (0.02) 


161 (22 to 345) 


Staphylococcus*** 


0.56 (<0.01) 


259 (168 to 353) 


Corynebacterium*** 


0.50 (<0.01) 


2,857 (1,440 to 4,722) 



a The negative and positive responders to MC were identified using indicator analysis, which also produced the indicator values. 

b *, significant change in prevalence only; **, significant change in absolute abundance (i.e., load) only; ***, significant change in both prevalence and absolute abundance. 



Diagnostics, Indianapolis, IN) and stored at — 80°C until analysis. In this 
analysis, we evaluated the enrollment and year 1 swabs from 77 control 
and 79 intervention arm participants, selected at random from among all 
married men who, together with their spouses, remained persistently HIV 
negative during the trial. The study was approved by four institutional 
review boards: the Science and Ethics Committee of the Uganda Virus 
Research Institute (Entebbe, Uganda), the HIV subcommittee of the Na- 
tional Council for Science and Technology (Kampala, Uganda), the Com- 
mittee for Human Research at Johns Hopkins University's Bloomberg 
School of Public Health (Baltimore, MD), and the Western Institutional 
Review Board (Olympia, WA). 

Sample processing. We processed samples from each participant in 
the same batch to control for interrun variation. For each sample, we lysed 
100 fil of eluted transport medium using enzyme-free chemical and me- 
chanical lysis. We purified the lysate using a Qiagen AllPrep DNA/RNA 
minikit (Qiagen, Valencia, CA) and performed DNA elution using 100 u.1 
of buffer EB. Additional methodological details can be found in Text SI in 
the supplemental material. 

Bacterial load quantification and 16S rRNA gene-based pyrose- 
quencing analysis. We quantified the bacterial load, measured as the bac- 
terial 16S rRNA gene copies per /id of coronal sulcus swab eluent, using a 
broad-coverage qPCR assay (28). We also generated bar-coded V3-V6 
amplicons using broad-coverage fusion PCR primers, which were pooled 
and sequenced on the Genome Sequencer FLX instrument (Roche Ap- 
plied Science, Branford, CT). Resultant pyrosequences were chimera 
checked (29), demultiplexed, and quality checked (30). We performed 
taxonomic classification using the Ribosomal Database Project Naive 
Bayesian Classifier (RDP release 10, update 28) (31). Detailed description 
of the bioinformatics analyses can be found in Text S 1 in the supplemental 
material. 



After stringent filtering, pyrosequencing yielded a total of 104,425 16S 
rRNA gene sequences for samples from control men at enrollment and 
90,560 at year 1; for the intervention group, there were 88,834 16S rRNA 
gene sequences at enrollment and 66,265 at year 1 . These sequences rep- 
resented 18 phyla, 31 classes, 49 orders, 121 families, and 306 genera at a 
>80% bootstrap confidence level after excluding taxonomic groups with 
only a single sequence detected from the full sample set. For Clostridiales 
and Clostridiales family XI, many sequences could not be further classified 
at a >80% bootstrap confidence level. These were included in the data set 
as unclassified Clostridiales and unclassified Clostridiales family XI, re- 
spectively. 

Bacterial load comparison. We expressed the bacterial load change in 
each participant over time as a log response ratio (LRR) using the follow- 
ing: ln[(bacterial load at year l)/(bacterial load at baseline)] (32). LRR 
quartiles and means for participants from each group were plotted in the 
R (version 2.13.1) software environment (33) and compared used a two- 
tailed t test with unequal variance at an a value of 0.05. 

16S rRNA gene-based microbiota comparative analysis. We ana- 
lyzed the coronal sulcus microbiota based on operational taxonomic unit 
(OTU), i.e., the unique bacterial groups detected at each taxonomic level. 
We converted the per-participant OTU data into four metrics: prevalence, 
relative abundance, absolute abundance, and log-transformed absolute 
abundance. 

We calculated each OTU's prevalence as (total number of participants 
positive for the OTU in group X)/(total number of participants in group 
X) and the relative abundance as (number of sequences assigned to the 
OTU in participant A)/(total number of sequences from participant A). 
We calculated absolute abundance using the formula (relative abundance 
of each OTU in participant A) X (bacterial load in participant A) and the 
log-transformed absolute abundance as ln(absolute abundance + 1). 
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For the 40 most common genera at enrollment, we compared the 
baseline prevalences and relative abundances between the study arms us- 
ing the chi-square test and two-tailed f test, respectively. We assessed the 
change in prevalence (i.e., APrevalence) at year 1 in each arm using a 
presence-absence data matrix and a two-tailed paired t test. All P values 
were adjusted for false discovery. The APrevalences between the 
circumcised and uncircumcised men were further compared based on 
the following: AAPrevalence = [(^Prevalence^^^,^) — 
(APrevalence control )]. 

We compared the change in overall microbiota composition visually 
based on family-level log-transformed absolute abundance data using 
nonmetric multidimensional scaling (nMDS) and Bray-Curtis distance 
(34-36). The resultant nMDS plots were annotated with centroids and 
95% confidence ellipses (34). We assessed the microbiota change over 
time for each study group using permutational multivariate analysis of 
variance (PerMANOVA) (34) based on the log-transformed absolute 
abundance data in Euclidean distance. 

We assessed the change in microbiota biodiversity in each group using 
two biodiversity metrics: diversity (D), calculated as D = Simpson's di- 
versity index, and evenness (£), calculated as E = D/S, where S is richness 
(37). Evenness reflects the dominance by many (i.e., high evenness) versus 
few (i.e., low evenness) OTUs, whereas richness is a measurement of the 
total number of unique OTUs present. We calculated A£ and AD for each 
individual and applied bootstrapping to generate random control- 
intervention pairs (i = 1,000). We estimated the mean E and D effect sizes 
(ES) as follows: mean £ ES = mean (AE^^,^ xi - A£ control K ) and 
meanD ES = mean (AD interventionXi 

ADcontroi n)> as weu as me accom- 
panying 95% confidence intervals. 

We identified indicator bacterial genera impacted significantly by MC 
with indicator species analysis using log-transformed data (38). For these 
indicator genera, we quantified the mean MC effect sizes and the 90% 
confidence intervals. Detailed description of the statistical analyses can be 
found in Text S 1 in the supplemental material. 

Literature review. We performed a literature review of the oxygen 
tolerance of the 40 most common genera in the uncircumcised-group 
microbiota. Bergey's Manual of Determinative Bacteriology (39) was used 
for Corynebacterium spp., Lactobacillus spp., Staphylococcus spp., and 
Streptococcus spp. For others, we performed a search of the MEDLINE 
database via the PubMed tool with a cutoff date of April 2012 using a 
combined term of the applicable genus name and "nov" to identify pub- 
lications defining new species within the genus. Detailed results from the 
literature review can be found in Table S2 in the supplemental material. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi: 10.11 28/mBio.00076- 1 3/-/DCSupplemental. 

Text SI, DOC file, 0.1 MB. 

Figure SI, PDF file, 1.1 MB. 

Table SI, DOCX file, 0.2 MB. 

Table S2, DOCX file, 0.1 MB. 
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